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Abstract 

We analyze the semileptonic B c — > B*£ + £~ decay in the frame work of the Stan- 
dard Model. We calculate the B c to B* transition form factors in QCD sum rules. 
Analytical expressions for the spectral densities and gluon condensates are presented. 
The branching ratio of the B c — ► B*£ + £~ decay is calculated, and it is obtained that 
this decay can be detectable at forthcoming LHC machines. 
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1 Introduction 



The double heavy B c mesons have been firstly discovered by the CDF Collaboration [1], 
with mass m*jg = (6.4 ± 0.39 ± 0.13) GeV. Recently CDF Collaboration [2] announced an 
accurate determination of the B c meson mass, vtib c = (6.2857 ±0.0053 ±0.0012) GeV . The 
study of the B c meson has received a great interest, due to its special properties: firstly, 
its lowest bound state is composed of two heavy (charm and beauty) quarks with open 
flavor; secondly, this meson attracts the interest of physicists for checking predictions of 
the pertubative QCD in the laboratory; and lastly, the weak decay channels of B c meson 
are reacher compared to the corresponding decay channels of B q (q = u, d, s) and can be 
divided into three classes: 

• b quark decaying weakly, with the c quark as spectator, e.g., B c — > J/ip£v^ 

• c quark decaying weakly, with the b quark as spectator, e.g., B c — > B s £pf, 

• The annihilation channels like B c — > £u e , B c — > D*~K°*. 

There are quite a few number of theoretical works studying various leptonic, semileptonic 
and hadronic exclusive decay channels of B c meson (for a review, see [3]). 

From experimental point of view, the study of weak semileptonic decays of B c meson 
is quite important for the determination of Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements, leptonic decay constant /b c , etc. 

Much more B c mesons and more detailed information about their decay properties are 
expected at the forthcoming LHC accelerator. In particular, this holds true for the dedicated 
detectors BTeV and LHCB which are specially designed for the analysis of B physics where 
one expects to see up to 10 10 B c mesons per year with a luminosity 10 34 cm~ 2 s~ 1 [4] 

The rare decays constitute on of the most important class of decays for a more precise 
determination of the parameters of the SM, as well as looking new physics beyond the SM 
[5-7], since the flavor changing neutral currents (FCNC) are absent in the SM at tree level 
and appear only at loop level, due to the running of virtual particle in the loop. 

Of the FCNC processes involving K and B mesons, main attention has been focused on 
B° — B° mixing, b — > s£ + £~ , b — > 57, s — > d£ + £~, s — > duu, etc, since top quark runs in the 
loop. 

Study of the FCNC decays in charm sector has not received enough attention. This can 
be explained by the fact that in the SM, D° — D° mixing [8,9], as well as FCNC decays 
[10-12], are expected to be very small. Moreover, long distance effects are quite huge, since 
the loop in charm decay involves light down quarks. 

In the present work we present a detailed analysis of the semileptonic B c — > B^£ + £~ 
decay in the framework of sum rules method. This mode is likely to be observed in forth- 
coming accelerator experiments. Note that the radiative decay B c — > B*j, which at quark 
level described by c — > u + 7 transition is investigated in [13, 14]. 

The plan of this work is as follows: In the following section QCD sum rules of the three- 
point correlators are considered and sum rules for the form factors that are responsible for 
the B c — > Bll + £~ decay are constructed. In section 3 we present our numerical results and 
conclusions. 
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2 Sum rules for transition form factors 



In the present section we derive sum rules for the form factors that controls the B c — > 
B* a £ + £~ decay This decay is described by the c — > u£ + £~ transition at quark level. The 
matrix element for the c — > u£ + £~ decay can be written in the following form: 



M = 



Ga 



C , 9 // (m c )M7 At (l - 7 5 )c^£ + Ci (m c )u7 M (l - 75)^7^7^ 
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- 2m c C 7 (ra c )m<7„ i ,— (1 + 75)^7/ 



(1) 



At /i = m c , Cg^ is given by [15] 

C 9 e// (m c ) = C 9 K) + ]T A, 
where 
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i7r (for s > 4z 2 ) 



where s = s/m 2 , Zi = rrii/m c . The logarithmic dependence on the internal quark mass 
mj in T(zj) in Eq. (2) cancels a similar term in the function Fi(xi) entering in Cg(mw), 
leaving no spurious divergences in the limit m 8 — > 0. It should be noted here that QCD 
corrections do not effect C10, i.e., Cio(m c ) = Cio(mw). 

The QCD corrections are particularly important for the Wilson coefficient C7 and in the 
numerical analysis we use the two-loop Q CD-corrected value of C7 which is calculated 
in [11]. 

The values of the Wilson coefficients at /i — mw are given by the following expressions 
(see for example [10, 13]) 

C 7 (m w ) = - ^i F 2{xi) , 

i=d,s,b 



1 



C 9 (m w ) = — Yl \ i [(C box (x i ) + C z (x i ))-2s 2 w (F 1 (x i ) + C 2 

% i=d,s,b 

C 10 (m w ) = ~ J2 K(C box (x t ) + C z ( Xl )) , 

S W i=d,s,b 



X, 



(3) 



where Aj = V ci V* i: Xi = mj/mly. The functions -Fi(x), F 2 (x), C box (xi) and C z (xi) are those 
derived in [16] and are all given in Appendix-A. 

Similar to the b — > s£ + £~ transition, the Wilson coefficient in the c — » u£ + £~ transition 
receives long distance contributions which have their origin in the real qq intermediate state, 
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i.e., p and u> mesons. These contributions can be written via the following replacement in 
C e 9 ff (m c ) [17] 



r eff , 3?r „ myTy^f+e- 

CT , my. — S — IVdyX Vf 



where my and Yy t are the resonance mass and width. The Fudge factor Ki is determined 
in [10] to have the values k p = 0.7 and k u = 3.1. 

Having the matrix element for c — > u£ + £~ transition at hand, our next problem is the 
calculation of the matrix element for the B c — > B*£ + £~ decay. It follows from Eq. (1) 
that, in order to calculate the amplitude of the B c — > B^£ + £~ decay, the following matrix 
elements are needed 

(S;|« 7m (1-7 5 )c|B c ) , 

These hadronic matrix elements of the B c — > B*£ + £~ decay can be parametrized in terms 
of the form factors in the following way: 



(s;0/,e) |« 7m (1 - 75 )c| B c {p)) = e^ a pe*»p a p 



m Bc + m B * 



- i ^(m Bc + m^JA^q 2 ) - (e*g)P M 



(e*q) 2 -^[A 3 (q 2 )-A (q 2 )]q, 



A 2 (q 2 ) 



m Bc + m B * 



(4) 



(S;(p', e) \ua^(l + 75 )c| S c (p)> = 2ie lipaf) e*>> P a P ' f, T 1 tf) + [e*(m| c - m|.) - (e*g)^]T 2 (g 2 ) 



9m 



m^ c — m 



2 



T 3 (g 2 ) , 



(5) 



where = (p + p')mj Q/j, — (p ~ p')h an d £* is the polarization 4-vector of the vector B* 
meson. Note that the form factor A 3 (q 2 ) can be written as a linear combination of Ai and 
A 2 as: 



2m B * 



2m B * 



^A 2 (q 2 ) , 



(6) 



and in order to guarantee the finiteness of the results at q 2 = 0, A 3 (0) = A o (0) should be 
satisfied. 

The identity 



o> 7 5 



a/3 



leads to the following relation between the form factors Ti and T 2 , 

T\(0)=T 2 (0). 
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Moreover, it is shown in [18] that the form factors Ti(q 2 ) and A (q 2 ) are also related, i.e., 



A (q 2 ) = -Uq 2 ) , 
and so are T 3 (g 2 ), A 1 (q 2 ) and A 2 (q 2 ) through the relation 



T 3 (q 2 ) = - 



2m B * 2m B *(m Bc + m B *) 



(7) 



It follows from these relations that in describing the B c — > B*£ + £~ decay, we need to know 
only V, A\ and A 2 . In further analysis the values of these form factors at q 2 = are needed 
(see section-3). Note that, at this point the form factor 21(0) is calculated in [14], and 
hence we do not present its expression in this work. 

As has already been noted, in order to calculate the these form factors appearing in the 
B c — > B*£ + £~ decay, we explore the three-point QCD sum rules [19]. 

For the evolution of the above-mentioned form factors in the frame work of the QCD 
sum rules, we start with the following three-point correlation functions: 

IW = -/ dl'xd i ye^^(0\T{j Biv (y)J lt {0)J Be {x)} |0> , (8) 

where J B » = b{y) r ) v u{y) is the interpolating current of the B* meson, J M = -u(0)7 M (l + 
75)c(0), and J Bc = c(x)i(l + 75)6(2;) is the interpolating current of the B c meson. The 
contribution of the scalar part of this current to B c is zero. 

The Lorentz structures of these correlators can be written in the following forms: 

n^ +A) = e vll afi p a p' P U v + Tl Al g^ + U 2 V^ V + U A2 V u q^ + IlgP^ + IL^ . (9) 

The phenomenological part of these correlators can be calculated by inserting a complete 
set of intermediate states with the same quantum numbers the currents J B * and J Bc posses 
in the correlation functions (4) and (5), and expressing these functions as the sum of the 
contributions of the lowest lying and excited states, we get 

KvJpIpI I 2 ) = ~ J ;\ Bc) ^ Jgc|0) + contributions from higher states .(10) 

The matrix elements in (10) are defined as follows: 

t {B c \-c{i +1 ,m = ^^. (11) 

Performing summation over the polarization of B* meson on the matrix elements (B*\ii'y tl (l+ 
•y 5 )c\B c ) that are given in Eq. (4), we get for the physical part: 

= fBjB*m 2 Bc r fi 2V 

(m b + m c ) (pi - m 2 Bc ) (p\ - m 2 B « ) \ vw * m Bc + m B * 
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+ i{m Bc + m B .) - + Mi 

\ B* J 

_^_, + ^) ( , 3 _, o) }. (12) 

The expressions of the form factors V, Ai and A 2 can be determined by choosing the 
coefficients of the Lorentz structures e u ^ a pp a p 113 , and V^qv in the correlator function 
U ( y +A \ respectively. 

On the other side, the three-point correlator function can be calculated by operator 
product expansion (OPE) in the deep Euclidean region p\ <C (mj + m c ) 2 , pi <C m 2 . 

The time ordered products of currents in the three-point correlator function in Eq. (1) 
can be expanded in terms of a series of local operators with increasing dimension, as is 
shown in the following: 

-J d 4 xd 4 ye-^ x ^T{.J Biu .^J Bc } = (C0UI+ (C 3 ) ull qq + (C^G^G^ 

+ (C 5 )^qa a ^G a/3 q + {C^qVqqV'q , (13) 

where {C^) Vfl are the Wilson coefficients, / is the unit operator, G a/3 is the gluon field 
strength tensor, T and V are the matrices appearing in the calculations. Considering the 
vacuum expectation value of the OPE, the correlator function can be written in terms of 
the local operators as: 

Tl^(pipiq 2 ) = (CoW + (CsMqq) + (C 4 ) Ufl (G 2 ) + (C 5 ) u ,(qa a(S G^q) 

+ (C 6 ) ulM (qTqqT'q) . (14) 

The values of the heavy quark condensates are related to the vacuum expectation value of 
the gluon operators in the following manner: 

<««-i4v (G!, -i4' (c3) ' (15) 

where Q is the heavy quark and the heavy quark condensate contributions are suppressed 
by inverse of the heavy quark mass, and for this reason we can safely omit them. 

It should be noted here that, in principle, the light quark condensate does give con- 
tribution to the correlator function, but its contribution becomes zero after double Borel 
transformation. Therefore, the only nonperturbative contribution to the above-mentioned 
correlator function comes from the gluon condensate. 

As a result, in the lowest order of perturbation theory, the three-point functions receive 
contribution from the bare-loop and gluon condensates, i.e., 

Tl.iplpW) = Ur(plpiq 2 ) + uf\plplq 2 )^(G 2 ) . (16) 

7T 

The bare-loop contribution can be obtained using the double dispersion representation 

1 r p per (s,s',Q 2 ) 
jjper _ _ — / — 1 — 55 d s ds' + subtraction terms , (17) 

(2ir) 2 J (s - p 2 )(s' - j/ 2 ) 
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in the variable p 2 and p' 2 , where Q 2 = —q 2 . The integration region in Eq. (17) is determined 
by the inequalities 



-1 < 



2ss' + (s + s' + Q 2 )(m 



mi 



s) + 2sm\ 



<+l , 



(18) 



A 1 /2( S;S ',_g2) A l/2( m 2 )m 2 ?s ) 

where A (a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 26c. 

The spectral densities p^ er (s, s', Q 2 ) can be calculated by using the Gutkovsky rule, i.e., 
by replacing the propagators with Dirac-delta functions 

1 



k 2 



m c 



-2in5(k 2 - m 2 ) . 



(19) 



After standard calculations, we get for the spectral densities: 
Pv 



2N c m c 



Pa 1 = 



A 3 / 2 (s,s',-Q 2 ) 
2N r m r 



A 3 / 2 (s,s',-Q 



(2s' Ai - uA 2 ) , (20) 
m c s' + m 2 m 2 (u — 2s') + m 2 s'{u — 2s) + Q 2 {m\u — ss' — ml) 



+ l -\( s ,s',-Q 2 ) 



PA 2 



N c m c 



\y 2 {s,s',-Q 2 ) 
where iV c = 3 is the color factor, u 



- (2sA 2 - uAi) +B 1 + C-D-E 



(21) 



(22) 



s + s' + Q 2 , Ai = s + mf - m 2 c , A 2 = s' + mf, and 



explicit forms of the functions B, C, D, and E are given in appendix-B. According to the 
QCD sum rule philosophy, contributions coming from the excited states are approximated 
as bare-loop contribution, starting from some thresholds s and s', in accordance with the 
quark-hadron duality. Note that we neglected 0(a s /7r) hard gluon corrections to the bare 
loop diagrams, since they are not available yet. However, we expect their contributions to 
be about ~ 10%, so that if the accuracy of QCD sum rules is taken into account, these 
corrections would not change the results drastically. 

The next problem is calculation of the gluon condensate contributions to the correlation 
functions. The gluon condensate contributions to the three-point sum rules are described 
by the diagrams presented in Fig. (1). The calculations of these diagrams are carried out 
in the Fock-Schwinger fixed-point gauge [20-22] 



, 



where A a ^ is the gluon field. In calculating the gluon condensate contributions, integrals of 
the following types are encountered: 



I [a,b,c] = j 
7 M [a,6,c] = J 
I^[a,b,c] = 



d 4 k 


1 


(2tt) 4 [jfe2 


— m 2 ] a [{p + k) 2 — m 2 ] 6 [{jpf + k) 2 ] c 


d 4 k 




(2tt) 4 [jfc2 


— m 2 ] a [(p + k) 2 — m 2 f [{jpf + k) 2 ] c 


d A k 


kpk v 



(2tt) 4 \k 2 - m 2 \ a [(P + k) 2 - m 2 c } b [{p> + kf 



(23) 
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These integrals can be calculated by continuing to Euclidean space-time and using 
Schwinger representation for the Euclidean propagator 



+ m A 1 (a) Jo 



(24) 



k 2 + m 2 r(a) 

which is very suitable for the Borel transformation since 

B ir 2(M 2 )e- ap2 = 5(1/ M 2 - a) . (25) 

Performing integration over loop momentum and over the two parameters which we have 
used in the exponential representation of propagators [21], and applying double Borel trans- 
formations over p 2 and p' 2 , we get the Borel transformed form of the integrals in Eq. (23) 
(see also [21]) 

/■I \a+b+c 

h(a, b, c) = , R J rf !, rfh , rfA {M^ 2 ~ a ' b (M 2 ) 2 - a - c V(o(a + b + c - 4, 1 - c - b) , 



i67r 2 r(a)r(6)r(c) 
i 

2 



4 (a, b, c) = ]- h(a, b, c) + I 2 (a, b, c) + - A (a, b, c) - I 2 (a, b, c) 



V(a, b, c) = J 6 (a, b, c)g^ + J (2/4 + 7 3 + A)?^ + \ ( - A + A)7>, 
+ -( - h + J 3 )P,g M + -( - 2/ 4 + h + h)q^v , 

where 

h(a,b,c 



(26) 



I 2 {a,b,c 
h(a,b,c 
Ii{a, b, c 
I 5 (a,b,c 
I 6 (a,b,c 



(_1 )a+b+c+l 

1- : \_! . {Ml) 2 - a -\Mlf- a - c U {a + b + c- 5,1 - c-b) , 



lQ7i 2 T(a)T(b)T(c 

_]\a+6+c+l 



i67r 2 r(a)r(6)r(c 



J 327r 2 T(a)T(b)T(c 



lQ7r 2 T(a)T(b)T(c 



i67r 2 r(a)r(6)r(c 



{M 2 f- a -\Ml) 2 - a - c U Q {a + b + c-5,l-c-b) , 
(M 2 ) 2 - a - b (M 2 ) 4 - a - c U (a + b + c- 6,1 - c-b) , 
{Mlf- a - b (Mlf- a - c U {a + b + c- 6,1- c-b) , 
(M 2 ) 4 - a - b (M 2 ) 2 - a - c U (a + b + c-6,l-c-b) , 
(Mlf- a - b {Mlf- a - c U (a + b + c- 6,2 -c-b) , (27) 



16tt 2 T(a)T(b)T(c) 

where Mf and M| are the Borel parameters in the s and s' channel, respectively, and the 
function Uo(a, (5) is defined as 



roo 

Wo (a, b) = / dy(y + M\ + M 2 2 )V cxp 
Jo 



--B - B lV 
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where 





2 

m c 
Ml 




B-! = 


m 


B = 


1 




M\ Ml 


B 1 = 


m\ 


MlMl 



m 2 b Ml + Ml(ml + m 2 c ) 



(28) 



Hat in Eqs. (26) and (27) means that they are double Borel transformed form of integrals. 
Performing double Borel transformations on the variables p 2 and p' 2 on the physical parts 
of the correlator functions and bare-loop diagrams and equating two representations of the 
correlator functions, we get the sum rules for the form factors V, A x and A 2 : 



V 



x 



(m b + m c )(m Bc + m B *) 2 /Af? m %. m* 
— _ e Be' ie a « 

VBjn%jB*mB* 
1 f s 'o f s ° 

^L ds 'L pv{s ' s '' Q2)e 



s/M* e -s'/Ml _ i 



247T 2 \ 7T 



a. 



G 2 )C, 



A 1 



(m b + m c 



e m% c /Ml e ml t /Ml 



X 



e-^e-'M-i 



x 



fB c m 2 B jBzm B - l (m Bc + m B 

{i£>'/> (s ' s ' Q2) 

(m b + m c )(m Bc + m B *) 

fB c m 2 B J B *m B *(m 2 Bc + 3m 2 B , + Q 2 , 



a 



247T 2 \ 7T 



S G 2 )C^ 



24vr 2 \n )C * 



where sq and s' are the continuum thresholds in pseudoscalar B c and B* channels, respec- 
tively, and the lower bound integration limit of s is given as 



sl = 



(ml - Q 2 - m 2 - s'){m 2 c s' + m 2 b Q 2 ) 
K 2 + Q 2 )(m 2 h - s') 



(29) 



Explicit expressions of C\ , Cf 1 and Cf 2 are all presented in Appendix-C. 

At the end of this section we present the dilepton invariant mass distribution for the 
B c — > B*£ + £~ decay. Using the parametrization of the B c — > B* transition in terms of form 
factors and Eq.(l), the matrix element of the B c — > B*£ + £~ decay can be written as 



M 



Got 



m Bc 



(30) 



where 



4 = G 1 (s)e, pa ^p a p' (i - tG 2 (s)e; + tG 3 (s)(e*q)V, + iG A {s){e* q)q^ , 

Jl = ffi(s)e^W - tH 2 (s)e; + iH 3 {s){e*q)V, + iif 4 (S)(e*g)g„ , (31) 



where = V^/m Bc , % = q f ,/m Bc and s = q 2 /m 2 Bc , and 



1 +f 
G 2 (s) = (1+r) 

GJs) 1 



2r, 



C e 9 ff A 1 (s) + ^(l-r)C e 7 ff T 2 (s) 



1 _p 



(l-f)C^ f A 2 (s) + 2f c C e 7 ff 



1-f 2 
T 3 (3) + — T 2 (s) 



G 4 (s) = -^ ff [(l + f)A 1 (s)-(l-f)A 2 (s)-2fA (s)} - 2f c C e 7 ff T 3 (s) 



:C w V(s) , 



H 2 (s) = (l+^doA^s) , 
H 3 (s) = ^C w A 2 {s) , 

H 4 (s) = lc l0 \(l + f)A 1 (s)-(l-f)A 2 (s)-2fA (s) 



(32) 



where f = m B */m Bc and f c = m c /m Bc . 

Using Eq. (30), the dilepton invariant mass distribution takes the following form: 



dT G 2 a 2 s 5 r- fl^A n ff\ ,2 l^ 2lrrl 



+ 
+ 



1 

4f 2 

4f 2 



A 



^! + 8f 2 (5 + 2 f 2 )) \G 2 \ 2 + (^(3 - v 2 ) + 8f 2 S^ |ff 



(3 - v 2 ) |G 3 | 2 + (hz - v 2 ) + 4f 2 (2 + 2f 2 - §)) \H 3 



2^2 j^ 3 - "'X 1 - f2 - *)R*[G 2 GS] + ^(3 - *; 2 )(1 -f 2 -s)+ 4f 2 A^Re[ff 2 iJ*] 



2A^f(Re[ff 2 ff 4 *] - (1 - f 2 )Re[# 3 # 4 *]) + ^sX |# 4 



(33) 



where 



A = A(l, f 2 , s), v 



N 



4m» 



is the lepton velocity, r 2 



3 Numerical analysis 

In this section we present our numerical calculation of the form factors Ai, A 2 , A , V, Ti, 
T 2 and T 3 . The values of the input parameters appearing in the sum rules for the form 
factors are: m Bc = 6.4 GeV, m B * = 5.325 GeV, f Bc = 0.385 GeV [24], f B * = 160 MeV 
[25], m b = 4.8 GeV, m c (fi = m c ) = 1-26 GeV, {(a s /Ti)G 2 ) = 0.012 GeV 4 [19]. The decay 
constants of B c and B* mesons are determined from the two-point QCD sum rules. As 
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has already been noted, in bare-loop calculations we neglect 0(a s /ir) corrections. For 
consistency, these corrections are also neglected in the calculations for the leptonic decay 
constants. 

The parameters s and s' Q , which are the continuum thresholds of B c and B* mesons, 
respectively, are also determined from the two-point QCD sum rules, and they are taken 
to be So = 50 GeV 2 and s' = 35 GeV 2 . These continuum thresholds are determined from 
the conditions that guarantees the sum rules to have the best stability in the allowed Mf 
and Mf region. 

The Borel parameters Mf and Mf in the sum rules are auxiliary parameters and physical 
quantities should be independent of them. Therefore it is necessary to look for a working 
regions of these parameters where physical results exhibit best stability. The working 
regions of Mf and Mf are determined by requiring that the continuum and higher states 
contributions are effectively suppressed, which ensures that the results do not sensitively 
depend on such excited states. We require also that the contribution gluon condensate 
is not too large, which guarantees that the contributions of higher dimensional operators 
are small. Our analysis verifies that the working regions 20 GeV 2 < Mf < 40 GeV 2 and 
10 GeV 2 < Mf < 15 GeV 2 of the Borel parameters satisfy both of the above-mentioned 
requirements. 

In these regions of Mf and Mf , the gluon condensate contribution constitutes approxi- 
mately 5%, and higher state contributions constitute at most 30% of the total result. 
Our numerical results for the form factors at q 2 = are: 

V(0) = 0.09 ±0.02 , 
Ai(0) = -0.17 ±0.03 , 
A 2 (0) = 1.10 ±0.20 , 

71(0) = T 2 (0) = -A)(0) = 0.23 ± 0.04 , (34) 

and the value of the form factor T 3 (0) can easily be obtained from Eq. (7). 

The errors are estimated by the variation of the Borel parameters Mf and Mf , the 
variation of the continuum thresholds So an d s' , the variation of b and c quark masses and 
leptonic decay constants f Bc and f B *. The main uncertainty comes from the thresholds 
and the decay constants, which is about ~ 20% of the central value, while the other un- 
certainties are small, constituting a few percent. Note that all the uncertainties are added 
quadratically. Here we would like to make the following cautionary note. It is well known 
that for heavy quarkonia, where the quark velocities are small, the a s /v corrections caused 
by the Coulomb-like interaction of quarks, are essential, where v is the quark velocity. In 
our case, we have two expansion parameters, a s /v 1 and a s /v 2 , where v 1 and v 2 are the 
relative velocities of quarks (be) and (bu) (for massless u quark v 2 = 1). When these cor- 
rections are taken into account the value of the form factors at Q 2 = are twice as greater. 
In further numerical analysis we will omit the dependence of the form factors on q 2 , which 
gives small contributions to the overall result. Indeed, the maximum value of q 2 in the 
decay under consideration is about (m Bc — m B *) 2 ~ 1 GeV 2 and assuming that the simple 
pole model correctly describes the q 2 dependence of the form factors, it is easy to see that 
q 2 /m 2 Bc is about 1/36; i.e., the results could be changed maximally about 3%. 

Integrating Eq. (32) over q 2 in the whole physical region and using the total mean life 
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time r ~ 0.46 ps of B c meson [26], the branching ratio of the B c — > B*£ + £ decay is 

f 1.3 x 1(T 9 , 



B(B C - fl^r) = 



3.9 x 10 



-7 



where the upper value corresponds to the case when only short distance contributions are 
taken into account, and the lower one corresponds when short and long distance contribu- 
tions due to the p and u resonances are taken into account. 

It follows from this result that, the B c — > B*£ + £~ decay with the above-presented width 
can be measurable at LHC 

In conclusion, we calculate the B c — > B* transition form factors V(q 2 ), Ai(q 2 ), A 2 (q 2 ) 
and Ti (q 2 ) in the framework of QCD sum rule. Our calculations show that this rare, semilep- 
tonic decay can be measurable at LHC. Using the results of the form factors calculated in 
the present work, we estimate the branching ratio of the B c — > B*£ + £~ decay. 
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Appendix-A 



In this appendix we present the expressions of the functions F ± , F 2 C box and C z , which 
enter to the expressions of the Wilson coefficients C 7 (m w ), C 9 (m w ) and C w (m w ): 



C box ( Xl ) 
C z { Xl ) 

Fifa) = Q d 
+ 



1 Xj In Xj 

+ 



Xi-i (xi - iy 



x^ 3 1 3 2x ^ x^ _ 
7 ~ %x~^l + 8 8(xi - l) 2 * ' 



1 1 13 

+ 



\2xi-l 12(^-1) 2 2(^-l) 3 
2 1 (2 1 5 1 



X i 



7 1 



+ 



1 1 



\3 (xi - 


I) 2 


6 (xi - 


1)3 


13 1 




1 1 




12 ( Xi - 


I) 2 


2( Xi - 


1)3J 


35 1 




5 1 





1 1 
+ 7:1 



X; 



InXi 



.L ; 



+ 



1 1 



Q Xl -l 12( Xl -iy 6(^-1)3 2(^-1) 



X{ 111 X{ , 



F 2 (x i ) = -Q d | 
'1 1 



1 3 

+ 



+ 



1 1 3 

4xi-l + 4(x,-l) 2 ' 2(.r, - 1)3 
9 1 3 1 



3 x 2 In 

2( Xl -iy 



3 

+ ~- 



2^-1 4(x,-l) 2 2( Xi -l) 



3 In Xi 

Xi ~~ " 



2 2(^-l) 4 



(A.l) 



where Xi = mf/m^, and Qd is the down quark charge. Note that , in these expressions, we 
omit the gauge dependent terms 7(£, a^) [16], because these terms are canceled out in the 
combinations C box ( Xi ) + C z ( Xi ) and Fi(xj) + C z ( Xi ). 
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Appendix-B 

In this appendix we present the expressions of the functions B 1 , C, D, and E, which 
appear in the calculations of the spectral density p A2 (see eq. (22)). 



A{S, S , — J 

+ 2m 2 s'[g 4 - (s - s') 2 + 3m c 2 (« - 2s') + Q 2 u 

+ s' 2 [6m 4 + Q 4 + (s - s') 2 + 6m 2 (« - 2s) + 2Q 2 (s - 2s')] } , 



C = D 



+ s' 



I {ml(u - s') [ - g 4 + m 2 b (2Q 2 - s) + s s 



+ s' 



+ 2m^ 



£ = 



A 2 M',-Q 2 y 

s(2Q 2 - s)(Q 2 + s) + m 4 (g 2 + 2s) - m 2 (Q 4 + 6Q 2 s + s 2 )" 
m 4 + m 2 (s - Q 2 ) - s(Q 2 + 2s)] + s' 3 (m 2 - s) - 3m 4 s'w 

- m 2 b ((u - s') 2 + s'(s - g 2 ) - 2s') - s'(Q 2 (u - 2s) - 2s 2 + s'u)] } , 

L——Us 2 [m 4 c (u + 6ss') + 2m 2 s[Q 4 - (s - s') 2 + Q 



A 3/2( S; s',-Q 2 Y 

+ m 4 [g 4 + (s - s') 2 + 2Q 2 (-2s + s')] + 2m 2 [ - ml(q 2 (u - 2s) - 2s 2 + s'u) 

+ s(-(u- s') 2 + s'(-s + g 2 ) + 2s' 2 )] } , (B.l) 



where 



4A 1 /2( S;S ',-g2) 
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Appendix-C 



In this appendix we give the explicit expressions of the coefficients of the gluon condensate 
which enter to the sum rules for the form factors V, A 1 and A 2 , respectively. 

C\ = 192m c 7 1 (l,3,l) + 192m^7i(l,4,l)-32m c / 1 (2,l,2) 

+ 64m c 7i(l, 1, 3) + 32m c 7i(2, 1, 2) + 128m c m 2 ii(2, 1,3)- 128m c 7j°' 1] (2, 1, 3) 
+ 64m c 7!(3, 1, 1) + 192m c m^7!(3, 1,2) - 128m c 7{°' 1] (3, 1,2) + 64m c m£ii(3, 1,3) 
- 128m c m^7T' 1] (3, 1, 3) + 64m c 7{ 0,2] (3, 1, 3) + 32m c A(l, 2, 2) + 32m c A(2, 2, 1) 
+ 32m c m 2 7 1 (2, 2, 2) - 32m c 7} 0,1] (2, 2, 2) - 32m c 7i(l, 2, 2) + 192m c m 2 ii(4, 1, 1) 
+ 96m c /i(2,2,l) , (C.l) 

= 96m c m 2 7 (l,3, 1) - 96m c 4°' 1] (l, 3, 1) + 96m^7 (l, 4, 1) 

- 96m;Lf{ ' 1] (l,4, 1) - 384m c 7 6 (l,3, 1) - 384m^7 6 (l, 4, 1) - 16m c 7 (l, 1, 2) 

- 16m c 7 (2, 1, 1) - 16m c m 2 7 (2, 1, 2) + 16m c 4° ' 1] (2, 1, 2) - 64m c 7 6 (2, 1, 2) 

- 16m c 7 (l, 1, 2) + 96m c m 2 7 (l, 1, 3) - 96m c i5 0,11 (l, 1, 3) - 16m c 7 (2, 1, 1) 

+ 48m c m 2 7 (2, 1, 2) - 80m c 7i°' 1] (2, 1, 2) + 96m c m^7 (2, 1, 3) - 1927n c 7n 2 b i [ ^ 1] (2, 1, 3) 
+ 96m c 4°' 2] (2, 1, 3) + 64m c m 2 7 (3, 1, 1) - 96m c 4°' 1] (3, 1, 1) + 64m c m^7 (3, 1, 2) 

- 160m c m27i°' 1] (3, 1, 2) + 96m c 7,f' 2l (3, 1, 2) + 32m c m^7 (3, 1, 3) - 96m c m^7i°' 1] (3, 1, 3) 
+ 96m c m 2 4° ,2] (3, 1, 3) - 32m c 4°' 3] (3, 1, 3) + 64m c 7 6 (2, 1, 2) - 128m c m 2 7 6 (2, 1, 3) 

+ 128m c 4°' 1] (2, 1, 3) - 256m c m 2 7 6 (3, 1, 2) + 128m c 4°' 1] (3, 1, 2) - 128m c m^7 6 (3, 1, 3) 
+ 256m c m 2 4°' 11 (3, 1, 3) - 128m c 4 ' 1] (3, 1, 3) + 32m c m 2 7 (l, 2, 2) - 32m c 7 ( f' 11 (i ! 2, 2) 

- 32m c 4°' 1] (2, 2, 1) + 16m c m^7 (2, 2, 2) - 32m c m 2 4° ,1] (2, 2, 2) + 16m c 7,f ,2] (2, 2, 2) 

- 128m c 7 6 (l, 2, 2) + 64m c m 2 7 6 (2, 2, 2) + 64m c 7f 1] (2, 2, 2) - 16m c 7 (l, 2, 1) 

- 16m c m 2 7 (l, 2, 2) + 16m c /{° 2, 2) - 64m c 7 6 (l, 2, 2) + 96m c m 2 7 (3, 1, 1) 

+ 96m c m£7 (4, 1, 1) - 96m c m^o 0,1] (4, 1, 1) - 384m c m 2 4(4, 1, 1) + 48m c 7 (l, 2, 1) 

+ 48m c m 2 7 (2,2,l)-48m c 7 ( f' 1] (2,2,l)-192m c 7 6 (2,2,l) , (C.2) 

= 96m c 7 2 (l, 3, 1) + 96m%{l, 4, 1) - 96m c 7 3 (l, 3, 1) - 96m%(l,A, 1) 

+ 96m c 7 5 (l, 3, 1) + 96m%(l, 4, 1) + 16m c 7 2 (2, 1, 2) - 16m c 7 3 (2, 1, 2) + 16m c 7 5 (2, 1, 2) 

- 32m c 7 (l, 1, 3) - 32m c 7 (2, 1, 2) - 32m c m 2 7 (2, 1, 3) + 32m c 7<f ' 1] (2, 1, 3) - 32m c 7 2 (l, 1, 3) 

- 16m c 7 2 (2, 1, 2) + 32m c 7 2 (3, 1,1) + 96m c m 2 7 2 (3, 1,2)- 64m c 7i°' 1] (3, 1, 2) 

+ 32m c m^7 2 (3, 1, 3) - 64m c m 2 4°' 1] (3, 1, 3) + 32m c 4°' 2] (3, 1, 3) + 16m c 7 3 (2, 1, 2) 

- 32m c m 2 7 3 (2, 1, 3) + 32m c 7 3 ' 1] (2, 1, 3) - 64m c m 2 7 3 (3, 1, 2) + 32m c 7 3 ' 1] (3, 1, 2) 
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- 32m c m 4 b I 3 (3, 1,3) + 64m c m 2 / 3 ' 11 (3, 1,3)- 32m c / 3 °' 2] (3, 1,3)- 16m c / 5 (2, 1, 2) 
+ 32m c m 2 / 5 (2, 1,3)- 32m c 4 0,11 (2, 1,3) + 64m c m 2 / 5 (3, 1,2)- 32m c 4 0,11 (3, 1, 2) 
+ 32m c m 4 b I 5 (3, 1, 3) - 64m c m 2 4°' 1] (3, 1, 3) + 32m c 4°' 2] (3, 1, 3) + 16m c J (l, 2, 2) 

- 32m c m 2 b I (2, 2, 2) + 48m c f 2 (l, 2, 2) + 16m c / 2 (2, 2, 1) - 48m c m 2 f 2 (2, 2, 2) 

- 16m c / 2 ' 1] (2, 2, 2) - 32m c / 3 (l, 2, 2) + 16m c mg/ 3 (2, 2, 2) + 16m c / 3 ' 1] (2, 2, 2) 
+ 32m c / 5 (l, 2, 2) - I6m c m 2 b l 5 (2, 2, 2) - 16m c 4°' 1] (2, 2, 2) - 32m c /i(l, 2, 2) 

+ 16m c / 2 (l, 2, 2) - 16m c / 3 (l, 2, 2) + 16m c / 5 (l, 2, 2) + 96m c m 2 / 2 (4, 1, 1) 

- 9Qm c m 2 b I 3 (A, 1,1) + 96m c m 2 / 5 (4, 1, 1) + 48m c / 2 (2, 2, 1) - 48m c / 3 (2, 2, 1) 
+ 48m c / 5 (2,2,l) , 

where 

(m 2 ) 1 (m 2 ) J I n (a,b,c) . 



a, 6, c) = (M 2 



M 2 2 



d" d J 
d(M?Yd(M2Y 
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